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Aerospace  Corporation  and  McDonnell  Douglas  Astronautics  Corporation. 
The  study  was  designed  to  develop  a  pogo  stability  model  yielding  eigen¬ 
values  and  eigenvectors  consistent  with  flight  data  over  a  range  of  vehicle 
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analyses  was  required.  Based  on  this  model,  it  was  possible  to  optimize 
the  accumulator  resistance  for  vehicle  stability. 


; 

■ 

1 


I 

I 

! 


UNCLASSIFIED 


4KCOMTV  CU  AMtFIC— TlOM  OF  THIS  Ml  CfttarMQ 


S.  Rubin,  Senior  Staff  Engineer 
Vehicle  Integrity  Subdivision 
Vehicle  Engineering  Division 


♦ 


I 


Publication  of  this  report  does  not  constitute  Air  Force  approval  of  the 
report ls  findings  or  conclusions.  It  is  published  only  for  the  exchange  and 
stimulation  of  ideas. 


L  £_ 

William  D.  Knox,  Colonel,  USAF 
Asst  Program  Director 
Expendable  Launch  Vehicles  SPO 

i 

i 


m 


i 

s. 

i 


CONTENTS 


I.  INTRODUCTION . .  .  .  .  5 

II.  CORRECTION  OF  MATHEMATICAL  MODEL .  13 

A.  Changes  to  Stability  Equations .  13 

B.  Changes  to  input  Data  . .  14 

1.  Analytical  Estimation  of  Lox  and  Fuel  Pump- 

Inlet  Compliance . 15 

2.  Inference  of  Fuel  Pump-Inlet  Compliance  from 

Flight  Data  . . 15 

III.  STABILITY  PREDICTIONS .  23 

IV.  TRANSFER  FUNCTIONS . 29 

V.  MODIFICATIONS  TO  LOX  ACCUMULATOR .  33 

A.  Optimization  of  Resistance .  33 

B.  Definition  of  Resistance  .  35 

VI.  SUMMARY  AND  CONCLUSIONS .  45 

REFERENCES  .  51 

APPENDIX . 53 

SYMBOLS .  61 


FIGURES  (Continued) 


18.  Effect  of  Superimposed  Random  Flow  on  Accumulator 

Resistance . 42 

19.  Maximum  Gimbal  Block  Accelerations  for  ITOS-F 

Relative  to  Other  Delta  Vehicles . 43 

20.  Comparison  of  Guidance -Compartment  Acceleration  for 

Five  ITOS  Vehicles .  49 

A-l.  Schematic  of  Delta  System  . . . . .  54 


SECTION  I 


INTRODUCTION 


Since  February  1972  The  Aerospace  Corporation  has  been  involved  in 
a  joint  study  with  the  McDonnell  Douglas  Astronautics  Corporation  (MDAC)  in 
an  effort  to  suppress  pogo  oscillations  on  a  new  class  of  Delta  vehicles,  in¬ 
volving  an  extended -long -tank  Delta  with  an  RS-27  engine.  The  culmination 
of  that  portion  of  the  study  reported  herein  has  been  a  successful  demonstra¬ 
tion  of  suppression  capability  on  a  prior  class  of  vehicle,  namely,  the  long- 
tank  Delta  with  an  MB -3  engine.  The  demonstration  was  achieved  on  the 
ITOS-F  vehicle,  and  the  means  of  suppression  was  a  resistive  lox  accumula¬ 
tor.  The  purpose  of  this  report  is  to  describe  the  studies  which  resulted  in 
the  successful  pogo  suppression  on  the  ITOS-F  flight. 

The  current  Delta  family  of  vehicles  was  preceded  by  the  Thor/Agena 
and  Thor/Delta  families.  Figure  1  depicts  this  greater  class  of  vehicles  in 
terms  ox  the  three  sizes  of  the  first  stage:  short -tank,  long -tank,  and  ex¬ 
tended-long-tank  configurations.  All  vehicle  configurations,  beginning  with 
the  inception  of  Thor/Agena  flights  in  late  I960,  have  exhibited  a  major  pogo 
instability  during  the  twenty  seconds  prior  to  first-stage  main -engine  cutoff 
(MECO).  Long-tank  and  extended-long -tank  Delta  vehicles  have  also  exhibited 
two  earlier  periods  of  pogo  activity. 

The  term  "pogon  is  a  nickname  given  co  self-excited  vibrations  exper¬ 
ienced  by  liquid  rocket  vehicles  arising  from  interactions  of  the  vehicle 
structure  with  the  propulsion  system  (Refs.  1  and  2);  various  studies  of  the 

^"Prevention  of  Coupled-Structure  -  Propulsion  Instability,"  NASA  Space 
Vehicle  Design  Criteria  (Structures),  NASA  SP-8055,  October  1970. 

^"Ftnal  Technical  Report,  Thor  Block  2  Test  Program  (20  cps  Oscillation 
Problem)."  R-3128,  North  American  Aviation,  Inc.,  August  1961. 
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phenomenon  are  reported  in  Refs.  3  through  8.  Figure  2  depicts  the  pogo 
events  for  selected  vehicles  from  each  class  of  Delta  vehicle  in  terms  of 
the  variation  of  longitudinal  acceleration  amplitude  at  the  first-stage  engine 
gimbal  block.  The  ITOS-F  results  should  be  compared  to  typical  results  for 
two  earlier  vehicles  in  its  class  (long-Unk  Delta  with  an  MB-3  engine):  (1) 
the  IDCSP/A  having  no  lox  accumulator,  and  (2)  the  TD-1  having  a  lox  accu¬ 
mulator  but  without  added  flow  resistance  at  the  entrance  to  the  accumulator 
as  on  the  ITOS-F.  The  reduction  in  the  severity  of  all  three  pogo  events, 
most  dramatically  the  MECO  event,  is  evident. 

From  mid-1963  to  mid- 1964,  a  panel  of  personnel  from  MDAC,  Rocket- 
dyne,  Lockheed  Missile  and  Space  Company,  and  The  Aerospace  Corporation 
("Thor  20-CPS  Panel")  provided  a  formal  mechanism  for  communication  re¬ 
garding  the  MECO  pogo  instability  on  the  short-tank  Thor/Agena  vehicles. 

3 

Radovcich,  N.A.,  "Analytical  Model  for  Missile  Axial  Oscillation  Induced  by 
Engine  -Struct ure  Coupling,  "  A658522,  Lockheed  Missiles  and  Space  Company, 
March  1965. 

4 

Davis,  W,  F.,  T.  F.  Lynch,  and  T.  R.  Murray,  "Thor  20  Cycle  Longitudinal 
Oscillation  Study,"  The  Shock  and  Vibration  Bulletin,  .34(2),  pp.  177-196, 
December  1964. 

^Woodward,  G.  B.,  and  D.  H.  Gunderson,  "Evaluation  of  High  Compliance 
Device  to  Suppress  20  cps  Longitudinal  Oscillations  in  the  Thor  Vehicle,  " 
TER-65-2,  North  American  Aviation,  Inc.,  July  1965. 

^Davis,  W.  F.,  et  al,  "Fe  sibility  Study  of  Compliant  Devices  to  Eliminate 
Longitudinal  Instability,"  SM-49005,  Douglas  Aircraft  Company,  August  1965. 

"Summary  Report,  Development  of  a  Longitudinal  Oscillation  Model  for  the 
Thor  MB-3  Propulsion  System,"  R-6337,  North  American  Aviation,  Inc., 
October  1965. 

g 

Lewis,  W.,  "Simplified  Analytical  Model  for  Use  in  Design  of  Pump-Inlet 
Accumulators  for  the  Prevention  of  Liquid-Rocket  Longitudinal  Oscillation 
(POGO),"  NASA  TN  D-5394,  1969. 


In  particular.  Aerospace  disseminated  information  obtained  from  the  Titan/ 
Gemini  vehicle  which  was  the  first  demonstration  of  successful  pogo  suppres¬ 
sion  (Ref.  1).  In  large  measure,  the  success  on  Titan/Gemini  was  attribu¬ 
table  to  a  ground  test  program  of  the  propulsion  system  which  led  to  its  satis¬ 
factory  mathematical  modeling  (Ref.  9).  Unfortunately,  no  such  test  program 
had  been  conducted  for  the  Thor  propulsion  system.  Over  the  years  extensive 
data  provided  statistical  measures  of  the  pogo  structural  loads.  These  loads 
were  accommodated  by  redesign  of  the  aft  structure  of  the  Agena,  by  requali¬ 
fication  of  Agena  equipments  to  higher  sinusoidal  vibration  levels,  and  by 
design  of  payloads  to  withstand  the  loads.  A  program  requirement  to  sup¬ 
press  the  instability  did  not  exist  and  so  studies  were  terminated.  Subsequent 
success  with  an  accumulator  on  the  Titan  IIIB  program  (Ref.  1)  prompted 
studies  of  similar  devices  for  the  short-tank  Thor/Agena  (Refs.  5  and  6).  It 
was  concluded  that  stability  could  be  attained.  Without  a  program  require¬ 
ment.  vehicle  modification  was  not  implemented. 

Beginning  in  1966,  long-tank  vehicles  replaced  the  short-tank  ones.  On 
long -tank  Thor /Delta  vehicles  two  earlier  pogo  events  were  first  observed 
(as  typified  by  the  IDCPA/A  vehicle  shown  on  Fig.  2).  Again,  statistics  were 
compiled  and  no  action  was  taken  to  correct  the  vehicle.  In  1970  MDAC, 
studying  the  pogo  stability  of  the  extended -long -tank  Delta  vehicle,  predicted 
higher  pogo  acceleration  levels  at  MECO  than  those  occurring  on  the  long- 
tank  vehicles.  It  was  evident  that  some  type  of  corrective  device  was  re¬ 
quired,  and  it  was  concluded  that  a  slightly  modified  version  of  the  device 
employed  on  the  S-II  Stage  of  the  Saturn  V  i  stalled  in  the  lox  suction  line 
would  stabilize  the  vehicle.  Aerospace,  requested  to  express  an  opinion, 
concurred  that  a  lox  accumulator  would  probably  result  m  a  stable  vehicle  at 
MECO  based  on  the  close  parallel  of  the  structure/propulsion  dynamic  charac 
teristics  (as  then  understood)  of  the  Delta  vehicle  to  those  of  the  Titan  IIIB 
vehicle  which  has  been  successfully  stabilized;  no  comment  could  be  made  for 
the  two  earlier  pogo  events  for  which  no  parallel  existed  on  Titan, 

9  "  ’  — ~ 

'Wagner,  R.G.,  "Titan  II  Engine  Transfer  Function  Test  Results,"  Report 
No.  TOR -0059(6471  )-9,  The  Aerospace  Corporation,  January  1971. 
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The  MDAC  recommendation  to  install  the  modified  Saturn  accumulator 
or.  extended -long -tank  Delta  vehicles  was  accepted.  In  order  to  demonstrate 
feasibility  of  the  concept,  an  accumulator  was  flown  on  a  series  of  long-tank 
vehicles.  The  flight  data  from  these  long-tank  vehicles  (as  typified  in  Fig.  2 
by  the  TD-I  results  relative  to  those  from  IDCSP/A)  showed  not  a  stable 
vehicle  as  predicted,  but  a  shift  to  a  later  time  of  the  MECO  instability  and 
the  continued  presence  of  the  so-called  pre-MECO  instabilities. 

After  several  flights  with  the  lox  accumulator.  Aerospace  was  funded 
by  NASA/Goddard  to  conduct  a  parallel  study  with  MDAC.  Thic  report  sum¬ 
marizes  the  Aerospace  effort  in  improving  the  understanding  of  pogo  stability 
of  Delta  vehicles  with  the  MS-3  engine.  This  improved  insight  led  to  a  modi¬ 
fication  of  the  lox  accumulator  which  was  responsible  for  reduced  oscillations 
observed  during  the  ITOS-F  flight.  The  studies  reported  herein  have  been 
carefully  coordinated  with  similar  efforts  at  MDAC.  Additional  studies  are 
underway  for  the  extended-long -tank  vehicle  with  the  new  RS-27  engine,  but 
these  are  too  preliminary  to  be  reported  at  this  time. 

In  the  absence  of  tests  to  determine  propulsion-system  characteristics, 
the  approach  has  been  to  rely  on  a  thorough  analysis  of  flight  data  in  an 
attempt  to  identify  the  defects  in  the  mathematical  model  (Ref.  1).  In  addi¬ 
tion  to  oscillograph  traces  to  spot  onset  and  cessation  of  periods  of  instability, 
power  and  cross -spectral  analyses  are  employed  to  yield  frequency  content  as 
well  as  relative  amplitude  and  phase  relationships  between  various  pressures 
and  accelerations.  The  periods  of  instability  are  related  to  the  eigenvalues 
from  closed-loop  stability  analyses,  a  positive  real  part  of  a  complex  eigen¬ 
value  indicating  an  unstable  system.  Likewise,  the  relative  amplitudes  and 
phases  among  various  flight-measured  parameters  are  directly  comparable 
to  the  eigenvectors  of  the  closed-loop  system.  Correction  of  the  stability 
model  is  based  on  improvement  of  such  correlations  by  changing  the  equations 
or  t.ie  input  data  or  both.  Section  II  discusses  the  changes  made  based  both 
on  flight  data  reduction  and  on  theoretical  modeling  refinements. 
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Sections  III  and  IV  discuss  the  changes  in  stability  predictions  and 
amplitude -phase  relationships  which  have  occurred  daring  the  study.  The 
result  was  a  good  understanding  of  the  MECO  pogo  event,  while  the  correla¬ 
tion  for  the  two  pre-MECO  events,  although  improved,  leaves  them  still 
clouded.  It  was  judged  possible  to  modify  the  accumulator  to  suppress  the 
MECO  instability,  without  deleterious  effect  on  the  pre-MECO  instabilities. 
The  ITOS-F  flight  displayed  success  with  a  substantial  reduction  in  the 
MECO  oscillations,  as  well  as  an  apparent  reduction  of  the  pre-MECO  events. 
Section  V  discusses  the  rationale  behind  the  accumulator  modification  and 
the  ground  testing  of  the  accumulator  required  to  validate  design  parameters. 
Section  VI  summarizes  the  effort  to  date  and  recommends  future  action. 


i 
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SECTION  II 


CORRECTION  OF  MATHEMATICAL  MODEL 


Initial  Aerospace  stability  analyses  of  the  Delta  vehicle  were  conducted 
using  a  set  of  stability  equations  developed  for  the  Titan  vehicle  as  described 
in  Refs.  10  and  11.  It  was  recognized  that  the  equations  were  not  completely 
applicable  to  the  Delta  vehicle  but  nonetheless  it  was  felt  that  only  minor 
modifications  would  result  in  a  model  which  would  describe  the  required 
phenomena.  As  the  study  progressed,  comparisons  of  stability  predictions 
with  flight  data  made  it  apparent  that  the  form  of  the  model  was  not  wholly 
adequate,  and  some  extensive  modifications  were  made.  Additionally,  the 
data  used  to  compute  the  coefficients  of  the  equations  were  continuously  up¬ 
dated  throughout  the  study.  This  process  of  changing  the  model  and  com¬ 
paring  it  to  the  flight  experience  was  carried  out  in  three  primary  stages. 

A.  CHANGES  TO  STABILITY  EQUATIONS 

The  first  stage  of  change  to  the  equations  involved  a  modification  to 
include  effects  of  oscillatory  outflow  from  the  propellant  tanks.  A  derivation 
of  the  equations  describing  tank  potential  outflow  and  an  investigation  of  the 
potential  significance  for  stability  had  been  undertaken  during  investigations 
of  the  Space  Shuttle  vehicle  (Ref.  12).  This  change  resulted  in  new  equations 
for  tank-bottom  pressure  and  in  changes  to  the  generalized  force  portion  of 
the  structural  modal  response  equation. 


Rubin,  S.,  "Instability  Model  of  Missile  Oscillation  Due  to  Propulsion  Feed¬ 
back,  "  Report  No.  TOR-269(4126)-28,  The  Aerospace  Corporation,  Sep¬ 
tember  1964. 

^Payne,  J.  G.,  "Pogo  Stability  Analysis  Program,"  Report  No.  TOR-0172 
(21221-10,  The  Aerospace  Corporation,  January  19?2. 

> 

'Rubin.  S. ,  R.  G.  Wagner,  and  J.  G.  Payne,  "Pogo  Suppression  on  Space 
Shuttle  -  Earlv  Studies,  "  NASA  CR-2210,  March  1973. 
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The  following  changes  were  made  during  a  second  stage  of  change:  (1) 
the  fuel  feedline  portion  was  revised  to  include  the  fluid/'structui a!  coupling 
owing  to  the  area  change  (below  the  lox  aft  bulkhead),  (2)  the  lox  accumulator 
equation  was  modified  to  include  a  vertical  inertance  term  (for  partial  gas 
fill).  (3)  the  influence  of  the  elevation  change  along  the  discharge  line  was 
included,  and  (4)  horizontal  sections  were  incorporated  in  the  oxidizer  and 
fuel  feedlines  ahead  of  the  pumps. 

Tn  a  third  round  of  change,  equations  were  modified  to  account  for  the 
influence  of  lateral  and  rocking  pump  motions  on  flow  in  the  feedlines  and  on 
generalized  forces.  Also  implemented  was  a  more  refined  set  of  combustion 
relationships  (i. e.,  individual  combustion  delays  of  the  two  propellants). 

App  .dix  A  contains  a  complete  description  of  the  current  Aerospace 
pogo  stability  model.  It  should  also  be  pointed  out  that  the  MDAC  stability 
model  has  also  been  upgraded  (Ref.  13).  The  following  are  among  the  changes 
made  at  the  suggestion  of  Aerospace: 

1.  Representation  of  the  lox  suction  line  by  discrete  rather  than 
modal  coordinates  to  account  properly  for  feedline  and  accumu¬ 
lator  resistances 

2.  Inclusion  of  discharge  elevation  change 

3.  Account  for  oscillatory  tank  outflow 

B.  CHANGES  TO  INPUT  DATA 

For  initial  stability  analyses,  a  complete  set  of  stability -model  para¬ 
meters  was  obtained  from  MDAC  and  used  without  modification.  As  the  study 
progressed,  all  the  parameters  were  examined  by  both  MDAC  and  Aerospace 
and  few  numbers  have  not  been  changed.  For  example,  almost  all  inertance 
and  resistance  values  have  been  recomputed.  The  major  changes  to  the  data 
for  the  MB-3  engine,  however,  were  for  both  the  lox  and  fuel  pump-inlet 
cavitation  compliances.  A  recently  developed  analytical  estimation  technique 

^Leuner,  T.  R.  and  M.  J.  Morgan,  "MB-3  Powered  Long  Tank  and  Extended 
Long  Tank  Delta  Pogo  Stability  Analysis,  "  A3-250-ABD3-73-TM-3, 
McEtonnell  Douglas  Astronautics  Company,  May  1973. 


was  employed  for  both  pomps  (Ref.  14).  It  was  also  possible  to  infer  fuel 
pump-inlet  compliance  from  flight  data.  These  two  developments  of  com¬ 
pliance  are  discussed  in  the  next  sections. 

Figure  3  shows  the  original  and  modified  values  of  compliance  versus 
pump  cavitation  index,  with  the  compliance  expressed  in  terms  of  a  "bubble 
frequency,  "  f^,  given  by 

fb  *  25  V7^ 

By  definition,  f^  is  the  frequency  of  the  system  composed  of  the  line  inertance, 
L  ,  and  the  pump  compliance,  C.  .  As  explained  in  Ref.  10,  the  bubble  fre- 
quency  provides  a  useful  measure  of  the  influence  of  pump  compliance  in 
the  hydraulic  system. 

1.  Analytical  Estimation  of  Lox  and  Fuel  Pump-Inlet  Compliance 

The  analytical  estimation  of  lox  and  fuel  compliances  was  performed  by 
the  "stay  time"  method  developed  in  Ref.  14.  The  method  is  based  on  em¬ 
pirical  correlation  studies  of  oxidizer  and  fuel  compliances  determined  from 
ground  tests  of  the  Rocketdyne  H-l,  F-l,  and  J-2  turbopumps.  These  tests 
employed  sinusoidal  perturbations  in  flow  over  a  range  of  frequencies.  The 
results  for  the  MB-3  pumps  are  shown  in  Fig.  3  by  the  curves  labeled 
"analytical.  "  Note  that  they  are  much  lower  than  the  original  values  em¬ 
ployed  by  MDAC  at  the  higher  cavitation  indices.  In  view  of  these  discrep¬ 
ancies,  it  was  decided  for  Aerospace  to  undertake  an  independent  review  of 
the  flight  data  from  the  MB-3  engine. 

2.  Inference  of  Fuel  Pump-Inlet  Compliance  from 

Flight  Data 

A  means  of  obtaining  pump-inlet  compliance  is  to  determine  one  or 
more  propulsion-system  resonant  frequencies  from  test  data  and  then  infer 
that  value  of  compliance  required  in  the  analytical  model  to  match  the  observed 
__ 

Ghahremani,  F.  G.  and  S.  Rubin,  "Empirical  Evaluation  of  Fump  Inlet 
Compliance,"  Report  No.  ATR-73(7257)-l ,  The  Aerospace  Corporation, 
August  1972. 
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Fig.  3.  Bubble  Frequencies  for  Long-Tank 
Delta  Vehicle 


frequency  or  frequencies.  Employed  is  an  analytical  model  including  the 
feedline  with  an  accumulator  if  present,  the  line  termination  impedance 
provided  by  the  engine,  and  the  unknown  pump-inlet  compliance.  For  the 
short  lox  line  where  fluid  compressibility  is  negligible,  and  when  no  accumu 
lator  is  present,  it  is  necessary  to  determine  only  the  first  propulsion- 
system  resonant  frequency  in  order  to  infer  the  pump-inlet  compliance. 

With  a  lox  accumulator  of  known  compliance  present,  both  the  first  and  the 
second  resonant  frequencies  are  necessary  to  deduce  the  pump-inlet  com¬ 
pliance.  For  the  longer  fuel  feedline,  where  organ-pipe  effects  are  signifi¬ 
cant,  the  first  and  second  resonant  frequencies  must  be  identified  in  order 
to  deduce  both  the  pump-inlet  compliance  and  the  effective  acoustic  speed  in 
the  line. 

It  has  generally  proven  very  difficult  to  identify  propuls  ion -system 
resonant  frequencies  from  flight  data  of  Delta  vehicles.  The  technique  em¬ 
ployed  (Ref.  1)  is  to  perform  power  spectral  density  (PSD)  analyses  of  pump 
inlet  pressure  at  a  series  of  flight  times  and  to  attempt  to  identify  resonant 
frequency  variation  with  time  to  correspond  with  known  variation  in  cavita¬ 
tion  index.  PSD's  of  accelerations  must  also  be  obtained  in  order  to  identify 
peaks  due  to  structural  resonances  so  as  not  to  confuse  them  with  peaks  due 
to  propulsion  resonances. 

As  reported  in  Ref.  7,  spectral  results  of  pressure  and  acceleration 
data  were  available  from  the  short-tank  Thor/Agena  Vehicle  344.  Concur¬ 
rent  with  the  study  reported  herein,  MDAC  performed  PSD  calculations  on 
pressure  and  acceleration  data  from  three  long-tank  Delta  vehicles,  OSO-H, 
HEOS-A2  and  TD-1,  as  well  as  on  data  from  Telesat-A,  an  extended -long - 
tank  Delta  vehicle.  Additionally,  MDAC  obtained  spectrogram  displays 
(Ref.  1)  sometimes  called  "voiceprints",  of  the  data  for  several  of  these 
vehicles. 

All  of  the  spectral  data  mentioned  were  examined  by  Aerospace. 
Attempts  to  detect  the  necessary  lox  resonant  frequency  information 
for  inference  of  pump-iniet  compliance  were  unsuccessful.  For  the  fuel 
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system,  only  the  PSD's  from  the  short -tank  Vehicle  344  were  judged  to  indi¬ 
cate  the  first  and  second  fuel  resonant  frequencies  versus  time  (Fig.  4).  it 
was  this  information  which  served  as  the  basis  for  our  inference  of  fuel 
pump-inlet  compliance.  Some  com\nents  should  be  made  regarding  the 
Aerospace  interpretation  of  this  Vehicle  344  data.  The  first  fuel  system 
resonant  frequency  was  inferred  f/  mi  peaks  in  pressure  not  correctable 
with  peaks  in  structural  acceleration.  On  the  other  hand,  the  second  reso¬ 
nant  frequency  is  based  on  sporadically  visible  peaks  in  the  fuel  suction 
pressure  spectra  which  are  almost  entirely  matched  by  peaks  in  correspond¬ 
ing  acceleration  spectra.  Our  interpretation  is  that  such  peaks  are  only 
visible  in  fuel  pressure  spectra  when  a  structural  resonance  happens  to  be 
coincident  such  that  the  fuel  mode  responds  to  a  greater  degree.  The  con¬ 
fidence  with  which  these  peaks  can  be  identified  varies  with  flight  time  and 
so  the  influence  of  the  time  varying  cavitation  index  is  used  to  aid  in  identi¬ 
fying  the  expected  trend  in  frequencies  with  time. 

Also  shown  in  Fig.  4  is  the  Rocketdyne  first  fuel  suction  system  fre¬ 
quency  versus  flight  time  as  given  in  Ref.  7.  It  is  our  conclusion  that  the 
curve  from  Ref.  7  tracks  the  first  structural  mode  and  was  misinterpreted 
as  the  fuel -system  frequency.  Additionally,  early  use  of  this  data  did  not 
account  for  the  influence  on  the  observed  frequency  of  the  impedance  of  the 
engine  termination. 

In  order  to  infer  the  compliance  and  fluid  compressibility,  it  was 
necessary  to  construct  a  detailed  model  of  the  fuel  suction  system  as  it 
existed  on  Vehicle  344.  Figure  5  shows  the  form  of  the  model  of  the  suc¬ 
tion  line  and  Table  I  lists  the  engine  parameters  used  in  this  detailed  model. 
Pump-inlet  compliance  and  fluid  bulk  modulus  were  varied  in  order  to  match 
the  two  fuel  frequencies.  The  second  frequency  is  highly  dependent  on  bulk 
modulus  and  a  modulus  value  (B)  of  0.43  x  10  Ib/in.  produced  the  best 
match.  An  implication  of  using  this  value  of  modulus,  which  is  well  below- 
published  values,  is  that  the  second  fuel  resonance  for  the  long-tank  vehicle 
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is  then  in  the  30-Hz  region  and  can  play  a  significant  role  in  pre-MECO 
stability.  MDAC  has  chosen  to  base  their  fuel  model  on  a  value  of  modulus 
within  the  range  cf  published  values,  resulting  in  a  second  fuel  resonance 
nearer  to  40  Hz. 

Table  I.  Engine  System  Parameters  Used  in  Inference  of 
Pump-Inlet  Compliance  and  Fluid  Compressi¬ 
bility  from  Vehicle  344 

2 

Combustion  resistance  =  0.9  sec/in. 

Pump  resistance  =  1. 17  sec/in.^ 

2 

Discharge  resistance  =  3.06  sec/in. 

2 

Discharge  inertance  =  0.  0458  sec/in. 

Combustion  time  lag  =  0.0027  sec 

Pump  gain  =  2.  0 


SECTION  III 


STABILITY  PREDICTIONS 

As  mentioned  previously,  initial  Aerospace  stability  predictions  used 
the  basic  Titan  stability  model  and  MDAC  data.  The  resulting  stability  pre¬ 
dictions  together  with  the  MDAC  results  current  at  that  time  are  presented 
in  Figs.  6  and  7.  It  can  be  seen  chat  the  initial  Aerospace  results  substan¬ 
tiated  the  MDAC  prediction  that  inclusion  of  the  lox  accumulator  would 
stabilize  the  vehicle  near  MECO.  At  pre-MECO  times,  the  two  models  are 
not  in  agreement  with  each  other  or  with  flight  experience.  The  flight  data 
indicates  two  pre-MECO  periods  of  instability  both  with  and  without  an 
accumulator  (Fig.  2).  The  MDAC  model  at  this  time  predicted  one  instability 
with  no  accumulator  and  one  period  of  minimum  margin  with  accumulator. 

The  Aerospace  model  predicted  no  instabilities,  two  periods  of  minimum 
margin  with  no  accumulator,  and  a  single  minimum-margin  period  with  an 
accumulator. 

Figures  8  and  9  show  current  stability  predictions  from  both  the  Aero¬ 
space  and  MDAC  models  (Ref.  13).  Both  models  now  show  excellent  agree¬ 
ment  with  flight  results  near  MECO  in  that  the  incorporation  of  a  lox  accu¬ 
mulator  does  not  stabilize  the  vehicle  but  only  shifts  the  instability  nearer  to 
MECO.  No  substantial  improvement  is  evident  at  pre-MECO  times  except 
that  the  frequency  of  the  period  of  minimum  margin  in  che  MDAC  model  with 
an  accumulator  present,  and  in  the  Aerospace  model  with  no  accumulator 
present,  is  closer  to  the  frequency  of  the  second  pre-MECO  instability  seen 
in  flight.  Still,  neither  model  predicts  two  instabilities,  or  even  two  periods 
of  minimum  margin,  for  both  a  vehicle  without  and  with  an  accumulator. 

A  natural  question  upon  comparing  Figs.  6  and  7  with  Figs.  8  and  9  is. 
"What  changes  to  either  the  stability  equations  and/or  the  input  data  have  led 
to  the  improved  correlation?"  Since  most  of  the  changes  were  made  in 
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Current  Aerospace  Stability  Predictions 


conjunction  with  other  changes  (i. e.,  the  equations  and  the  data  were  fre¬ 
quently  changed  together),  it  is  difficult  to  isolate  the  influence  of  individual 
changes.  At  MECO,  the  major  influence  is  the  change  to  the  feedline 
resonant  frequencies  resulting  from  changes  to  pump- inlet  compliances  and 
to  distributed  compliance  on  the  longer  fuel  line.  The  resulting  second 
lox  line  frequency  of  the  configuration  with  an  accumulator  is  now  coincident 
with  the  structural  frequency  near  MECO,  whereas  it  was  previously  not 
coincident.  Any  major  influence  of  the  fuel  system  on  vehicle  stability  near 
MECO  has  been  removed  by  a  substantial  lowering  of  the  first  fuel  frequency. 

Although  the  changes  to  the  stability  picture  near  MECO  are  explained 
by  the  changes  to  the  propuls  ion -system  frequencies,  little  can  be  said  about 
the  changes  earlier  in  time.  Subsequent  to  the  initial  analyses,  the  primary 
period  of  interest  was  MECO;  no  pre-MECO  analyses  were  conducted  fol¬ 
lowing  many  of  the  changes.  The  changes  to  the  propulsion -system  frequen¬ 
cies  may  have  a  significant  influence  on  pre-MECO  results  as  well  as 
MECO  results. 
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SECTION  IV 


TRANSFER  FUNCTIONS 


A  method  of  evaluating  accuracy  of  pogo  stability  models  is  comparison 
of  the  complex  mode  shapes  (eigenvectors)  obtained  from  the  closed-loop 
stability  analysis  with  amplitude  ratio  and  phase  relationships  obtained  from 
flight  data  (Ref.  J ).  As  pointed  out  in  the  reference,  one  of  the  principle 
advantages  of  closed-loop  stability  analyses  vis-a-vis  open-loop  analyses  is 
the  ability  to  use  the  resultant  mode  shapes  as  a  basis  of  evaluating  the  sta¬ 
bility  model.  Based  on  experience  with  the  Titan  vehicle,  it  was  felt  that  this 
technique  would  be  a  useful  tool  for  assessment  of  the  Delta  stability  model. 

Application  of  the  technique  requires  careful  and  extensive  flight  data 
reduction.  Power -spectral  densities  (PSD),  cross -spectral  densities  (XSD), 
coherence  functions,  and  transfer  functions  must  be  determined  for  para¬ 
meters  of  interest  at  a  sequence  of  flight  times.  The  transfer  functions  are 
the  desired  amplitude  and  phase  relationships;  our  experience  has  led  to  the 
criterion  that  these  relationships  should  be  considered  invalid  when  the  co¬ 
herence  is  less  than  0.8.  This  process  can,  in  general,  be  simplified  for 
evaluation  of  the  structural  modes  to  only  PSD  plots  of  structural  accelera¬ 
tions.  This  simplification  is  possible  because  the  structural  modes  are  very 
lightly  damped,  the  coherence  among  the  accelerations  is  essentially  unity, 
and  the  relative  phasing  (0  or  180  deg)  is  usually  known. 

For  the  Delta  vehicle,  it  has  generally  been  possible  to  obtain  valid 
relationships  (coherence  >  0.  8)  only  during  periods  of  instability.  In  con¬ 
trast,  for  the  Titan  vehicle,  it  has  generally  been  possible  to  determine  the 
ratios  of  both  suction  and  chamber  pressures  to  gimbal  block  acceleration  in 
the  first  longitudinal  mode  throughout  most  of  Stage  I  burn  as  shown  in 
Fig.  10.  It  is  uncertain  whether  the  character  of  the  system  is  inherently 
different  or  whether  different  data  reduction  techniques  are  responsible. 
Nonetheless,  the  transfer  function  information  obtained  for  the  Delta  vehicle 
has  been  quite  useful  in  evaluating  the  stability  model. 
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Fig.  10.  Suction  ant!  Chamber  Pressure  per  Unit  Engine 
Acceleration  in  the  First  Longitudinal  Mode  for 
the  Titan  III  Vehicle 
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The  earliest  model  evaluation  via  transfer  functions  revealed  that  the 
MDAC  and  Aerospace  models  were  in  poor  agreement  with  each  other  as 
well  as  with  flight  data.  Both  Aerospace  and  MDAC  made  a  series  of  changes 
to  their  stability  models  (the  Aerospace  changes  are  described  in  Section  VI). 
Another  change  was  to  correct  the  lox  suction  pressure  computed  in  the  model 
to  account  for  the  actual  location  of  the  flight  measurement  (above  the  lox 
elbow).  Following  the  changes,  the  two  models  yielded  results  in  good 
agreement  with  each  other  at  MECO. 

Relationships  of  flight  data  to  the  Aerospace  model  at  MECO  are  shown 
in  Figs.  11  and  12.  The  comparison  with  flight  data  is  somewhat  clouded 
because  of  the  sparseness  of  good  flight  data;  e.g.,  it  is  difficult  to  identify 
consistent  trends  versus  frequency  over  a  wide  frequency  range  as  was  done 
for  the  Titan  data.  Additionally,  obtaining  high  coherence  data  only  during 
periods  of  instability  means  that  nonlinear  effects  may  be  present.  It  is  of 
interest  to  note  that  the  data  from  the  Telesat  flight,  which  had  the  most 
careful  data  reduction  (contiguous,  more  narrow  time  slices),  produced 
t.^nds  in  all  the  transfer  functions  which  agree  with  those  from  the  stability 
models  (see  especially  Fig.  12).  The  correlation  is  considered  to  be  mod¬ 
erately  good  and  provides  a  basis  for  a  reasonable  level  of  confidence  in  the 
stability  model  for  the  MECO  period. 

At  pre-MECO  times  the  stability  model  produces  transfer  functions 
which  are  in  very  poor  agreement  with  the  flight  data.  It  must  be  recalled, 
however,  that  the  stability  predictions  during  these  periods  are  also  not  in 
good  agreement;  perhaps  whatever  model  changes  improve  the  stability  por¬ 
trait  will  improve  transfer  functions. 
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SECTION  V 

MODIFICATIONS  TO  LOX  ACCUMULATOR 


A.  OPTIMIZATION  OF  RESISTANCE 

Historically,  the  function  of  an  accrmulator  has  been  to  detune  the 
propulsion  system  from  the  structural  system.  Prior  to  the  first  flight  of 
a  Delta  vehicle  with  an  accumulator,  it  was  felt  that  the  particular  accumu¬ 
lator  would  accomplish  this  detuning;  an  approach  successfully  demonstrated 
on  the  Titan  and  Saturn  vehicles.  Unfortunately,  this  assumption  was  based 
on  an  inadequate  model  of  the  lox  suction  system. 


The  lox  feed  system  model  used  at  the  time  did  not  reveal  that  incorpo¬ 
ration  of  the  accumulator,  while  detuning  the  first  lox  mode  from  the  struc¬ 
ture.  would  introduce  a  second  lox  mode  that  would  lead  to  an  instability. 

Had  it  been  possible  to  locate  the  accumulator  nearer  to  the  pump  inlet,  this 
situation  w'ould  not  have  occurred.  Since  it  was  judged  impractical  to  re¬ 
locate  the  accumulator  nearer  fo  the  pump,  an  alternate  solution  was  re¬ 
quired.  The  approach  recommended  by  Aerospace  was  to  add  resistance  via 
orifices  to  the  entrance  of  the  existing  accumulator.  It  was  hoped  that  the 
resulting  damping  increase  to  the  second  feedline  mode  would  be  sufficient  to 
eliminate  the  instability. 

A  study  was  undertaken  to  define  an  optimum  accumulator  resistance 
by  studying  its  effect  on  the  feedline  flow  impedance  looking  toward  the  tank 
from  a  position  just  upstream  of  the  pump  cavitation.  This  impedance,  Z^, 
can  be  written  in  real  and  imaginary  terms 


Z^  =  R  +  iX 


(1) 


A  maximum  value  of  the  real  part,  R,  was  the  criterion  for  an  optimum 
accumulator  resistance.  Referring  to  Fig.  13,  this  line  impedance  can  be 
written 


Z»  =  Z~  + 
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Fig.  13.  Schematic  of  Lox  Propulsion  <5,  stem 


where  Z,  =  R-.  +  ift-'L,,  Z.  =  R.  +  itul,. ,  and  Z  =  R  +  ia>L  .  The  accumu- 
33  311  1  a  a  a 

lator  impedance  Z ^  reflects  the  fact  that  the  accumulator  compliance,  C 
is  negligible  :n  the  frequency  range  of  interest. 

For  convenience,  the  search  for  an  optimum  resistance  was  conducted 
at  a  constant  frequency  of  150  rad/sec  (23.9  Hz),  a  value  approximately  equal 
to  the  line  resonant  frequency  at  MECO  both  with  and  without  an  accumulator. 
A  plot  of  the  imaginary  part,  X,  versus  the  real  part,  R,  of  the  line  impe¬ 
dance  as  a  function  of  the  accumulator  resistance,  R  ,  is  shown  below  in 

a 

Fig.  14.  Curves  are  shown  for  two  values  of  the  accumulator  incrtance, 

2 

L  ;  the  value  L  -  0.  0006  sec/in.“  pertains  to  the  unmodified  accumulator, 
a  a 

li  can  be  seen  that  an  optimum  accumulator  resistance  does  exist  (i.e.,  it 

Z 

produces  a  maximum  of  R ),  a  value  of  R  =■  0.  32  sec/in.  ljeing  optimum 

2  2  a 

for  L  =  0.  0006  sec  /in.  .  The  effect  of  accumulator  mertance,  L  ,  can 

Si  it 

also  be  seen;  a  doubling  of  L  significantly  lowers  the  maximum  .alue  of 

K  but  produces  little  changes  in  the  optimum  value  of  I'.  .  Since  a  maximum 

<1 

of  R  is  desired,  any  increase  in  ir.ertance  over  the  value  for  the  unmodified 
n<  (  uo-.ilafor  should  be  kept  to  a  minimum. 


Fig.  14.  Effect  of  Accumulator  Resistance  on  Lox 
Feedline  Impedance 


That  the  optimum  as  defined  by  the  feedline  impedance  would  also 

be  an  optimum  in  terms  of  vehicle  stability  was  verified  by  performing 

stability  analyses  versus  flight  time  for  selected  resistance  values.  Since 

the  modifications  studied  were  expected  to  increase  accumulator  inertance 

2  2 

somewhat,  L  was  chosen  as  0.  001  sec  /in.  for  this  study.  The  results 
a 

are  shown  in  Figs.  15  and, 16.  These  figures  show  that  near  optimum  resis¬ 
tance  will  produce  a  major  improvement  in  stability  margin,  while  there  is 
no  value  of  added  resistance  which  will  produce  a  gain  greater  than  one  of 
the  two  configurations  previously  flown  (i.e.,  no  accumulator  or  unmodified 
accumulator).  Further,  the  optimum  obtained  via  stability  analyses  agrees 
very  closely  with  that  obtained  via  the  feedline  impedance  procedure. 

B.  DEFINITION  OF  RESISTANCE 

The  results  of  the  previous  section  show  that  it  is  possible  to  achieve 
an  optimum  vehicle  stability  condition  if  the  accumulator  can  be  modified  so 
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as  to  produce  a  certain  resistance.  Unfortunately,  the  prediction  of  resis¬ 
tance  due  to  purely  oscillatory  flow  as  a  function  of  fluid  properties  and 
flow-passage  geometry  is  not  a  well-defined  procedure. 

In  order  to*  minimize  inertance,  it  was  proposed  to  modify  the  accu¬ 
mulator  by  placing  a  thin  orifice  plate  containing  a  series  of  small  holes 
within  the  entrance  region  to  the  accumulator.  The  final  configuration  con¬ 
taining  this  so-called  "hatband"  is  shown  in  Fig,  17.  The  accumulator  modi¬ 
fication  for  the  ITOS-E  flight  did  not  have  the  benefit  of  a  ground  testing 
program  to  confirm  that  the  desired  resistance  was  achieved.  This  decision 
was  based  on  scheduling  constraints  and  was  rationalized  with  the  fact  that 
(1)  previous  ITOS  missions  had  flown  both  with  and  without  an  accumulator 
and  had  experienced  about  the  same  pogo  levels,  and  (2)  that  any  added 
accumulator  resistance  is  helpful  in  theory  (Fig.  15). 

Two  sources  of  data  for  estimating  resistance  were  known  to  Aerospace 
prior  to  TTOS-E  flight.  The  first  source  was  frequency-response  tests  of  the 
fuel  bladder -type  accumulator  now  flying  on  the  Titan  vehicle  (Ref.  15),  and 
the  second  source  was  acoustic  tests  above  1  kHz  directed  toward  suppres¬ 
sing  combustion  instability  (Ref.  16).  Both  tests  involved  oscillatory  flow' 
through  various  patterns  of  holes  similar  to  the  conditions  of  the  proposed 
modified  accumulator. 

The  two  sources  of  data  seemed  to  yield  contradictory  results.  In 
interpreting  the  Titan  data,  it  *>-is  assumed  that  a  major  portion  of  the  total 
resistance  was  due  to  orifice  sistance  as  opposed  to  bladder  deformation 
or  rubbing.  It  was  further  assumed  that  the  RA  product,  where  A  is  the  total 
open  area  of  the  pattern  of  holes,  should  be  a  constant.  With  this  interpre¬ 
tation  of  the  Titan  data,  and  allowing  for  scatter  and  uncertainty  in  the  data, 
a  range  of  0.  1  <  R^A  <  0.  32  was  derived. 

__ 

Zehnle,  R.,  "T-lII-M  Phase  1A  Frequency  Response  Test  -•  Fuel  Accumu¬ 
lator,"  TM  0472/20-2-68,  Martin  Marietta  Corporation,  June  1968. 

1  6 

Garrison,  G.  D.,  et  al,  "Suppression  of  Combustion  Oscillations  with 
Mechanical  Damping  Devices  -  Interim  Report,  "  PWA  FR  3299,  Pratt  and 
Whitney  Aircraft,  August  1969. 


The  Titan  data,  involving  liquid  (Aerozine)  flow  oscillations  in  the 
frequency  range  of  3  to  30  Hz,  showed  resistance  to  be  reasonably  indepen¬ 
dent  of  flow  amplitude.  However,  the  acoustic  data,  based  on  N-,  gas  flow  in  the 
2  to  8  kHz  range,  showed  both  a  linear  and  nonlinear  region.  The  linear 
acoustic  data  correlated  well  with  Rayleigh  drag  theory,  but  when  this  linear 
data  was  used  to  estimate  the  Titan  test  resistance,  an  orifice  resistance  an 
order  of  magnitude  lower  than  measured  on  the  Titan  tests  was  computed. 

Since  it  was  felt  at  the  time  that  orifice  effects  should  have  been  dominant 
for  the  Titan  configuration,  it  was  decided  to  base  the  design  on  the  experi¬ 
mental  Titan  data.  Examination  of  Fig.  16  reveals  that  the  choice  of  an 

•> 

open  area  of  0.  5  in.~  for  the  hole  pattern,  resulting  in  a  range  of  R  of  from 

0.2  to  0.  64  (corresponding  to  0.  1  <  R  A  <  0.  32),  will  be  optimum  for  stability. 

The  configuration  for  ITOS-E  contained  a  pattern  of  1 /16-in.  diameter  holes 

2 

in  a  0.  030-in.  hatband  to  yield  the  desired  0.5-in.  area. 

This  configuration  was  flown  on  the  ITOS-E  vehicle  and  resulted  in  a 

stability  profile  very  similar  to  that  of  a  vehicle  with  no  accumulator,  indi- 

2 

eating  that  the  resistance  obtained  with  the  0.5-in.  area  was  effectively 
sufficient  to  block  the  accumulator.  Concurrent  with  the  flight,  an  additional 
source  of  resistance  data  was  pointed  out  to  Aerospace  bv  NASA/Goddard 
personnel  (Refs.  17  and  18).  The  new  data  appeared  to  be  superior  tc  the 
previous  data  since  it  was  the  result  of  carefully  controlled,  well -documented 
experiments  with  two  liquids  (allowing  the  effect  of  viscosity  to  be  assessed), 
and  it  did  cover  the  frequency  range  of  interest  for  the  Delta  vehicle.  Unlike 
the  Titan  data,  there  were  no  other  potential  sources  of  losses,  and  there  was 
little  scatter  in  the  data. 


Thurston.  G.  B.,  and  C.  E.  Martin,  "Periodic  Fluid  Flow  Through  Circular 
Orifices,"  J.  A.  S.  A.,  25(1),  January  1953,  p.  26. 

* 

Thurston,  G.  B.,  "Periodic  Fluid  Flow  Through  Circular  Tubes,"  J.  A.  S.  A.. 
24(6).  November  1952,  p.  653. 
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The  data  of  Refs.  17  and  18  correlate  well  with  the  previously  evaluated 
acoustic  data  in  that  both  linear  and  nonlinear  regions  are  present.  The 
linear  portions  correlate  well  with  Rayleigh  theory  as  did  the  acoustic  data. 
Additionally,  it  was  discovered  that  when  the  resistance  per  unit  area  was 
plotted  versus  particle  velocity,  the  new  data  correlate  reasonably  well  with 
the  acoustic  data  in  the  nonlinear  region.  These  data,  when  applied  to  the 
ITOS-E  accumulator,  indicate  operation  in  the  nonlinear  region  for  flow 
amplitudes  corresponding  to  very  small  vehicle  acceleration  amplitudes; 
thus,  the  accumulator  would  be  effectively  blocked  during  flight.  In  retro¬ 
spect,  it  was  concluded  that  the  orifices  were  not,  in  fact,  the  major  source 
of  losses  for  the  Titan  accumulator  as  had  been  previously  supposed. 

Based  on  the  latest  data,  it  was  necessary  to  define  a  new  hatband 
configuration  for  optimum  stability.  It  was  decided  that  a  pulsing  test  of 
several  candidate  configurations  would  be  essential;  there  were  still  some 
questions  regarding  resistance;  e.g.,  the  data  in  Refs.  17  and  18  were  col¬ 
lected  with  single  orifices,  rather  than  orifice  a* rays,  and  without  the  super¬ 
imposed  random  flow  oscillations  which  characterize  the  flight  situation. 
Pulsing  tests  were  conducted  by  MDAC  with  Aerospace  providing  a  basis  for 
the  test  plan.  It  was  suggested  by  Aerospace  that  testing  include  not  only 
variations  in  sinusoidal  input  levels  (designed  to  evaluate  the  nonlinear  char¬ 
acteristics  displayed  in  Refs.  17  and  18)  but  also  determine  the  influence  of 
random  inputs  superimposed  on  the  sinusoids.  Hatband  configurations  having 
various  open  areas  and  several  liquids  (to  assess  the  influence  of  viscosity) 
were  used  in  the  tests. 

The  test  data  for  purely  sinusoidal  flows  correlated  well  with  that  of 
Refs.  17  and  18  at  frequencies  up  to  20  Hz  but  showed  an  unexplained  higher 
resistance  at  higher  frequencies.  In  general,  the  addition  of  random  inputs 
had  the  effect  of  providing  a  resistance  "floor,  "  as  depicted  in  Fig.  18.  One 
unexpected  trend  was  a  slight  decrease  in  resistance  with  increasing  vis¬ 
cosity;  the  lowest  viscosity  fluid  tested  (a  Freon)  had  a  viscosity  not  very 
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different  from  that  of  lox,  and  the  dependence  was  not  a  strong  one.  In 
summary,  it  was  felt  that  the  resistance  was  well  defined  by  the  tests. 

The  tests  showed  that  an  open  area  of  1.  5  in.^  would  produce  a  resis¬ 
tance  very  near  the  optimum  value  and  fairly  independent  of  flow  amplitude 
(Fig.  18)  up  ic  that  corresponding  to  a  gimbal-block  acceleration  amplitude 
of  about  1  g.  A  hatband  having  this  area  was  incorporated  into  the  ITOS-F 
vehicle  flown  on  6  November  1973.  Examination  of  oscillograph  traces  of 
filtered  gimbal-block  acceleration  from  ITOS-F  shows  clearly  that  the  levels 
during  the  MECO  and  second  pre-MECO  pogo  events  are  dramatically  reduced 
relative  to  levels  observed  on  previous  long -tank  vehicles,  and  the  level 
during  the  first  pre-MECO  period  is  unaffected  (Fig.  19).  Figure  20  shows 

the  envelopes  of  guidance -compartment  accelerations  near  MECO  from  five 

2 

1TOS  flights  and  shows  that  the  1.5-in.  resistive  accumulator  produced 

levels  lower  than  vehicles  with  or  without  an  accumulator.  (Recall  from 

2 

Section  V-B  that  the  0.5-in.  resistive  accumulator  produced  a  resistance 
sufficiently  high  to  effectively  block  the  accumulator.  ) 

It  cannot  be  concluded  with  certainty  whether  the  ITOS-F  vehicle  was 
stable  at  MECO  as  predicted  analytically.  However,  the  character  of  the 
data  does  not  indicate  the  classic  build-up  and  decay  generally  associated  with 
an  instability  but  appears  similar  to  low-level  random  response.  This  change 
in  character  was  also  observed,  but  to  a  lesser  degree,  during  the  second 
pre-MECO  period. 


Fig.  18.  Effect  of  Superimposed  Random  Flow  on 
Accumulator  Resistance 
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SECTION  VI 


SUMMARY  AND  CONCLUSIONS 

The  reductions  in  pogo  acceleration  levels  on  the  ITOS-F  Vehicle  were 
rendered  possible  by  improvements  in  the  quality  of  the  mathematical  sta¬ 
bility  model,  which  set  the  stage  for  the  concept  of  a  resistive  lox  accumu¬ 
lator.  The  improvement  in  the  mathematical  stability  model  has  resulted 
from  joint  MDAC  and  Aerospace  studies.  The  major  Aerospace  contribu¬ 
tions  are  enumerated  below. 

1.  The  correction  of  the  lox  and  fuel  pump  cavitation  compliances. 
Detailed  review  of  Thor/Agena  Vehicle  344  data  (Ref.  7)  suggested 
that  the  fuel  compliance  was  probably  an  order  of  magnitude 
higher  than  that  previously  provided  by  the  engine  contractor 
(Fig.  3).  This  was  supported  by  subsequent  estimation  cf  the 

fuel  compliance  using  the  empirical  method  developed  from  data 
for  the  Saturn  F-l,  H-l,  and  J-2  pumps  (Ref.  14).  From  a  sta¬ 
bility  standpoint,  the  result  was  that  the  fuel  contribution  to  the 
MECO  pogo  instability  changed  from  significant  to  minor;  con¬ 
sequently,  MDAC  halted  consideration  of  a  fuel  accumulator.  The 
empirical  estimation  method  also  provided  new  lox-pump  cavi¬ 
tation  compliances  which  were  not  substantially  different  from 
previously  accepted  values  at  low  NPSH,  but  provided  substan¬ 
tially  greater  values  at  high  NPSH.  This  lox  compliance  change 
i.s  therefore  potentially  significant  for  the  early  pogo  events. 

2.  The  detailed  mathematical  modeling  of  the  lox  and  fuel  feedlines 
by  extension  of  the  Titan  feedline  models  (Ref.  11).  The  exten¬ 
sions  accounted  for  the  presence  of  elbows  in  both  feedline''  and 
the  presence  of  the  area  change  of  the  fuel  feedline  after  it  exits 
from  the  lox  tank  bottom.  Moreover,  the  restudy  of  Vehicle  344 
data  also  suggested  that  the  fuel  organ-pipe  frequency  was  lower 
than  previously  assumed.  A  fuel-line  modeling  study  showed  that 
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the  fuel  bulk  modulus  would  have  to  be  below  the  range  of  pub¬ 
lished  values  to  account  for  the  flight  observations.  Use  of  the 
lower  fuel  bulk  modulus  was  shown  to  be  of  significance  for  the 
first  pogo  event.  This  is  believed  to  account  for  differences 
between  Aerospace  and  MDAC  stability  results  for  this  event. 

The  acceptance  by  MDAC  of  the  Aerospace  discrete-parameter 
lox  feedline  model,  in  lieu  of  their  representation  with  damped 
modes  of  the  suction  system,  was  essential  to  the  proper  con¬ 
sideration  of  the  resistive  accumulator. 

3.  The  utilization  of  spectral  analyses  of  flight  data.  In  accord  with 
Ref.  1,  Aerospace  recommended  and  MDAC  performed  cross - 
spectral  analyses  of  pressure  and  acceleration  data  from  as  many 
flights  as  possible.  Comparison  of  the  resulting  amplitude  and 
phase  relationship  among  the  pressure  and  acceleration  variables 
have  been  correlated  with  relationships  from  eigenvectors  ob¬ 
tained  from  the  stability  analyses.  This  has  led  to  new  estimates 
of  pump  dynamic  gains  and  the  conclusion  that  the  current  mathe¬ 
matical  model  is  adequate  for  the  MECO  stability  event,  but  still 
inadequate  for  the  two  earlier  events. 


The  modeling  of  the  structural  dynamics  of  the  high-pressure  discharge 
ducts  by  MDAC  offers  hope  of  explaining  the  pre-MECO  pogo  correlation. 

The  measured  resonant  frequencies  of  the  RS-27  discharge  ducts  during 
ground  vibration  tests  appear  to  correlate  with  the  instability  frequencies. 

The  flight  instrumentation  of  these  ducts,  together  with  spectral  reduction  of 
the  data,  can  be  vital  to  an  understanding  of  duct  behavior  and  its  role  in  the 
stability  picture. 


Understanding  of  the  pogo  stability  situation  on  Delta  vehicles  has  been 
greatly  hampered  by  the  lack  of  special-purpose  dynamic  testing  of  the 
propulsion  system  such  as  have  been  available  for  other  vehicles  (Ref.  1). 

In  particular,  it  has  only  been  possible  to  obtain  correctable  data  during  the 
pogo  events,  and  so  the  information  is  limited  to  the  associated  relatively 
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discrete  frequencies  for  those  events.  It  is  unclear  why  the  success  achieved 
by  Aerospace  in  analyses  during  the  Titan  vehicle  Stage  I  flight  in  obtaining 
correlatable  data  fairly  continuously  over  the  entire  flight  (Fig.  10)  has  not 
been  achieved  for  Delta  flights.  It  may  be  significant  that  the  Telesat  data 
reduction,  which  used  narrower,  contiguous  time  slices,  gave  trends  in  all 
the  eigenvector  relationships  consistent  with  stability  analysis  for  the  MECO 
event.  Studies  aimed  at  improving  the  data-reduction  procedures  should  be 
carried  out  during  analysis  of  data  from  future  vehicles. 
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The  corrected  mathematical  model  showed  that  the  MECO  instability 
with  an  accumulator  present  resulted  from  a  coincidence  of  the  second  lox 
feedline  mode  with  the  first  longitudinal  structural  mode.  The  frequency  of 
this  second  lox  mode  was  very  much  dependent  on  the  positioning  of  the 
accumulator  relative  to  the  pump,  a  distance  equal  to  about  one -third  of  the 
total  feedline  length.  A  significant  raising  of  the  second  lox  frequency  by 
relocating  the  accumulator  nearer  to  the  pump  was  considered  to  be  imprac¬ 
tical  because  of  the  extensiveness  of  the  hardware  impact.  The  alternative 
proposed  by  Aerospace  was  to  add  resistance  to  the  entry  into  the  accumu¬ 
lator  to  increase  the  damping  in  the  second  lox  mode.  A  preliminary  cri¬ 
terion  for  optimization  of  the  resistance  was  that  the  real  part  of  the  feedline 
flow  impedance  be  maximized  at  the  frequency  of  the  second  feedline  mode. 

It  was  also  shown  that  the  inertance  of  orifices  used  to  produce  the  resistance 
should  be  minimized;  that  is,  a  large  number  of  small  diameter  holes  is 
preferable  to  fewer  larger  holes.  Subsequent  stability  analyses  of  the  vehicle 
confirmed  that  the  greatest  stability  improvement  was  indeed  achieved  for 
that  value  of  resistance  optimized  via  the  feedline  impedance  approach. 


1 


The  decision  to  implement  an  optimized  resistive  accumulator  on  the 
ITOS-E  vehicle  did  not  permit  the  recommended  ground  testing  to  investigate 
orifice  resistance.  Reliance  was  placed  on  resistance  data  from  the  testing 
of  the  Titan  fuel  accumulator  (Ref.  15).  The  interpretation  of  this  data  was 
subsequently  shown  to  be  erroneous  and  resulted  in  an  overly  high  resistance 
as  demonstrated  by  the  ITOS-E  flight  results  (Fig.  2).  The  high  resistance 


-47- 


produced  an  effective  blocking  of  the  accumulator  and  the  ITOS-E  pogo  events 
were  undistinguishable  from  those  on  a  vehicle  without  an  accumulator.  Data 
reported  in  the  literature,  which  came  to  light  after  the  ITOS-E  accumulator 
was  configured,  showed  (1)  that  the  majority  of  the  Titan  accumulator  resis¬ 
tance  must  have  been  from  something  other  than  the  presence  of  orifices, 
most  probably  from  material  dissipation  within  the  rubber  bladders  or  blad¬ 
der  rubbing,  and  (2)  that  the  orifice  resistance  was  nonlinear,  with  the  value 
of  resistance  being  prr  portional  to  the  amplitude  of  flow'  oscillations.  It  was 
clear  that  sinusoidal  tesiing  was  mandatory  to  investigate  fluid  viscosity 
effects,  possible  interaction  among  neighboring  orifices,  and  the  presence 
of  random  flow  fluctuations.  The  basis  for  the  test  plan  was  developed  by 
Aerospace  and  the  tests  were  conducted  by  MDAC.  There  results  the  ITOS-F 

configuration  of  the  resistance  accumulator,  having  a  total  orifice  area  of 
2  2 

1.5  in.  as  opposed  to  the  0.5  in.  for  the  ITOS-E  configuration.  The  success 
of  pogo  suppression  on  the  ITOS-F  is  displayed  in  Figs.  2  and  20. 

After  many  years  of  major  pogo  instabilities  on  the  Thor  family  of 
vehicles,  a  sufficient  understanding  has  apparently  been  achieved  resulting 
in  an  ability  to  suppress  instabilities  near  MECO.  Such  understanding  is  still 
lacking  for  the  earlier  pogo  events  which  suggests  that  it  would  be  premature 
to  discontinue  the  studies. 
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Fig.  20.  Comparison  of  Guidance -Compartment  Acceleration  for 
Five  ITOS  Vehicles 
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APPENDIX 


CURRENT  AEROSPACE  POGO  STABILITY  EQUATIONS 


Figure  A-l  is  a  schematic  representation  of  the  Delta  system  followed 
by  a  current  set  of  Aerospace  stability  equations.  It  should  be  pointed  out 
that  the  subscript  "o"  is  associated  with  fuel  suction  system  variables  and 
the  subscript  “f"  with  lox  system  variables.  On  the  Titan  vehicle,  the  fuel 
line  is  the  short  line  represented  incompressibly  while  the  much  longer 
oxidizer  line  is  treated  as  compressible.  The  reverse  is  true  for  the  Delta 
vehicle,  i.e.,  the  lox  line  is  the  shorter.  However,  the  programs  for  the 
Titan  had  been  written  with  variables  subscripted  "o"  referring  to  that  line 
represented  as  compressible,  and  this  notation  was  maintained  rather  than 
reversing  all  the  "o"  and  "f"  subscripts.  References  10  to  12  provide  back¬ 
ground  information  on  the  bases  for  the  equations. 


FUEL  SYSTEM 


Fig.  A-l.  Schematic  of  Delta  System 
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Flow  Continuity  Across  Oxidizer  Pump 


Incompressible  Oxidizer  Suction  Line  from  Accumulator 
Location  to  Pump 
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Oxidizer -Pump  Performance 
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Incompressible  Fuel  Suction  Line  from  Accumulator 
Location  to  Pump 
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11.  Fuel-Pump  Performance 
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12.  Oxidizer -Accumulator  Performance 
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3  3.  Fuel -Accumulator  Perlcrmance 
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Flow  in  Fuel  Suction  Line  from  Tank  Bottom  to 
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It,.  Flow  in  Oxidizer  Discharge  Line 
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SYMBOLS 


a»P>  y.6  constants  in  combustion  equation 

£  modal  damping  ratio  (ND) 
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SYMBOLS  (Continued) 


weight  density  (lb/in.^) 
modal  displacement  shape  (ND) 
circular  frequency  (rad/sec) 

horizontal  portion  of  section  of  line  between  accumulator  and 
pump 

vertical  portion  of  section  of  line  between  accumulator  and 
pump 

accumulator 

bubble 

chamber 

discharge 

engine 

lox 

lox  line  flange 

fuel  line  flange 

refers  to  the  i**1  mode 

fuel 

pump 

suction  line,  suction  (pump  inlet)  conditions 
tank 

tank  bottom 

section  of  suction  line  between  tank  and  accumulator 
section  of  suction  line  between  accumulator  and  pump 


